\(\int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx\) [536]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (verified)
   Maple [A] (verified)
   Fricas [B] (verification not implemented)
   Sympy [F(-1)]
   Maxima [F]
   Giac [B] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 23, antiderivative size = 453 \[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\frac {2 (c d-b e) \sqrt {d+e x}}{c^2}+\frac {2 (d+e x)^{3/2}}{3 c}+\frac {\sqrt {2} \left (b^3 e^2-b^2 e \left (2 c d+\sqrt {b^2-4 a c} e\right )+c \left (a \sqrt {b^2-4 a c} e^2-c d \left (\sqrt {b^2-4 a c} d-4 a e\right )\right )+b c \left (c d^2+e \left (2 \sqrt {b^2-4 a c} d-3 a e\right )\right )\right ) \text {arctanh}\left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {2 c d-\left (b-\sqrt {b^2-4 a c}\right ) e}}\right )}{c^{5/2} \sqrt {b^2-4 a c} \sqrt {2 c d-\left (b-\sqrt {b^2-4 a c}\right ) e}}-\frac {\sqrt {2} \left (b^3 e^2-b^2 e \left (2 c d-\sqrt {b^2-4 a c} e\right )+b c \left (c d^2-e \left (2 \sqrt {b^2-4 a c} d+3 a e\right )\right )-c \left (a \sqrt {b^2-4 a c} e^2-c d \left (\sqrt {b^2-4 a c} d+4 a e\right )\right )\right ) \text {arctanh}\left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {2 c d-\left (b+\sqrt {b^2-4 a c}\right ) e}}\right )}{c^{5/2} \sqrt {b^2-4 a c} \sqrt {2 c d-\left (b+\sqrt {b^2-4 a c}\right ) e}} \]

[Out]

2/3*(e*x+d)^(3/2)/c+2*(-b*e+c*d)*(e*x+d)^(1/2)/c^2+arctanh(2^(1/2)*c^(1/2)*(e*x+d)^(1/2)/(2*c*d-e*(b-(-4*a*c+b
^2)^(1/2)))^(1/2))*2^(1/2)*(b^3*e^2-b^2*e*(2*c*d+e*(-4*a*c+b^2)^(1/2))+c*(a*e^2*(-4*a*c+b^2)^(1/2)-c*d*(-4*a*e
+d*(-4*a*c+b^2)^(1/2)))+b*c*(c*d^2+e*(-3*a*e+2*d*(-4*a*c+b^2)^(1/2))))/c^(5/2)/(-4*a*c+b^2)^(1/2)/(2*c*d-e*(b-
(-4*a*c+b^2)^(1/2)))^(1/2)-arctanh(2^(1/2)*c^(1/2)*(e*x+d)^(1/2)/(2*c*d-e*(b+(-4*a*c+b^2)^(1/2)))^(1/2))*2^(1/
2)*(b^3*e^2-b^2*e*(2*c*d-e*(-4*a*c+b^2)^(1/2))-c*(a*e^2*(-4*a*c+b^2)^(1/2)-c*d*(4*a*e+d*(-4*a*c+b^2)^(1/2)))+b
*c*(c*d^2-e*(3*a*e+2*d*(-4*a*c+b^2)^(1/2))))/c^(5/2)/(-4*a*c+b^2)^(1/2)/(2*c*d-e*(b+(-4*a*c+b^2)^(1/2)))^(1/2)

Rubi [A] (verified)

Time = 2.90 (sec) , antiderivative size = 453, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.174, Rules used = {838, 840, 1180, 214} \[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\frac {\sqrt {2} \left (b c \left (e \left (2 d \sqrt {b^2-4 a c}-3 a e\right )+c d^2\right )+c \left (a e^2 \sqrt {b^2-4 a c}-c d \left (d \sqrt {b^2-4 a c}-4 a e\right )\right )-b^2 e \left (e \sqrt {b^2-4 a c}+2 c d\right )+b^3 e^2\right ) \text {arctanh}\left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {2 c d-e \left (b-\sqrt {b^2-4 a c}\right )}}\right )}{c^{5/2} \sqrt {b^2-4 a c} \sqrt {2 c d-e \left (b-\sqrt {b^2-4 a c}\right )}}-\frac {\sqrt {2} \left (b c \left (c d^2-e \left (2 d \sqrt {b^2-4 a c}+3 a e\right )\right )-c \left (a e^2 \sqrt {b^2-4 a c}-c d \left (d \sqrt {b^2-4 a c}+4 a e\right )\right )-b^2 e \left (2 c d-e \sqrt {b^2-4 a c}\right )+b^3 e^2\right ) \text {arctanh}\left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {2 c d-e \left (\sqrt {b^2-4 a c}+b\right )}}\right )}{c^{5/2} \sqrt {b^2-4 a c} \sqrt {2 c d-e \left (\sqrt {b^2-4 a c}+b\right )}}+\frac {2 \sqrt {d+e x} (c d-b e)}{c^2}+\frac {2 (d+e x)^{3/2}}{3 c} \]

[In]

Int[(x*(d + e*x)^(3/2))/(a + b*x + c*x^2),x]

[Out]

(2*(c*d - b*e)*Sqrt[d + e*x])/c^2 + (2*(d + e*x)^(3/2))/(3*c) + (Sqrt[2]*(b^3*e^2 - b^2*e*(2*c*d + Sqrt[b^2 -
4*a*c]*e) + c*(a*Sqrt[b^2 - 4*a*c]*e^2 - c*d*(Sqrt[b^2 - 4*a*c]*d - 4*a*e)) + b*c*(c*d^2 + e*(2*Sqrt[b^2 - 4*a
*c]*d - 3*a*e)))*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[d + e*x])/Sqrt[2*c*d - (b - Sqrt[b^2 - 4*a*c])*e]])/(c^(5/2)*Sq
rt[b^2 - 4*a*c]*Sqrt[2*c*d - (b - Sqrt[b^2 - 4*a*c])*e]) - (Sqrt[2]*(b^3*e^2 - b^2*e*(2*c*d - Sqrt[b^2 - 4*a*c
]*e) + b*c*(c*d^2 - e*(2*Sqrt[b^2 - 4*a*c]*d + 3*a*e)) - c*(a*Sqrt[b^2 - 4*a*c]*e^2 - c*d*(Sqrt[b^2 - 4*a*c]*d
 + 4*a*e)))*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[d + e*x])/Sqrt[2*c*d - (b + Sqrt[b^2 - 4*a*c])*e]])/(c^(5/2)*Sqrt[b^
2 - 4*a*c]*Sqrt[2*c*d - (b + Sqrt[b^2 - 4*a*c])*e])

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 838

Int[(((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_)))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[g*
((d + e*x)^m/(c*m)), x] + Dist[1/c, Int[(d + e*x)^(m - 1)*(Simp[c*d*f - a*e*g + (g*c*d - b*e*g + c*e*f)*x, x]/
(a + b*x + c*x^2)), x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*
e^2, 0] && FractionQ[m] && GtQ[m, 0]

Rule 840

Int[((f_.) + (g_.)*(x_))/(Sqrt[(d_.) + (e_.)*(x_)]*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)), x_Symbol] :> Dist[2,
Subst[Int[(e*f - d*g + g*x^2)/(c*d^2 - b*d*e + a*e^2 - (2*c*d - b*e)*x^2 + c*x^4), x], x, Sqrt[d + e*x]], x] /
; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0]

Rule 1180

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
 q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rubi steps \begin{align*} \text {integral}& = \frac {2 (d+e x)^{3/2}}{3 c}+\frac {\int \frac {\sqrt {d+e x} (-a e+(c d-b e) x)}{a+b x+c x^2} \, dx}{c} \\ & = \frac {2 (c d-b e) \sqrt {d+e x}}{c^2}+\frac {2 (d+e x)^{3/2}}{3 c}+\frac {\int \frac {-a e (2 c d-b e)+\left (c^2 d^2+b^2 e^2-c e (2 b d+a e)\right ) x}{\sqrt {d+e x} \left (a+b x+c x^2\right )} \, dx}{c^2} \\ & = \frac {2 (c d-b e) \sqrt {d+e x}}{c^2}+\frac {2 (d+e x)^{3/2}}{3 c}+\frac {2 \text {Subst}\left (\int \frac {-a e^2 (2 c d-b e)-d \left (c^2 d^2+b^2 e^2-c e (2 b d+a e)\right )+\left (c^2 d^2+b^2 e^2-c e (2 b d+a e)\right ) x^2}{c d^2-b d e+a e^2+(-2 c d+b e) x^2+c x^4} \, dx,x,\sqrt {d+e x}\right )}{c^2} \\ & = \frac {2 (c d-b e) \sqrt {d+e x}}{c^2}+\frac {2 (d+e x)^{3/2}}{3 c}-\frac {\left (b^3 e^2-b^2 e \left (2 c d+\sqrt {b^2-4 a c} e\right )+c \left (a \sqrt {b^2-4 a c} e^2-c d \left (\sqrt {b^2-4 a c} d-4 a e\right )\right )+b c \left (c d^2+e \left (2 \sqrt {b^2-4 a c} d-3 a e\right )\right )\right ) \text {Subst}\left (\int \frac {1}{-\frac {1}{2} \sqrt {b^2-4 a c} e+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{c^2 \sqrt {b^2-4 a c}}+\frac {\left (b^3 e^2-b^2 e \left (2 c d-\sqrt {b^2-4 a c} e\right )+b c \left (c d^2-e \left (2 \sqrt {b^2-4 a c} d+3 a e\right )\right )-c \left (a \sqrt {b^2-4 a c} e^2-c d \left (\sqrt {b^2-4 a c} d+4 a e\right )\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {1}{2} \sqrt {b^2-4 a c} e+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{c^2 \sqrt {b^2-4 a c}} \\ & = \frac {2 (c d-b e) \sqrt {d+e x}}{c^2}+\frac {2 (d+e x)^{3/2}}{3 c}+\frac {\sqrt {2} \left (b^3 e^2-b^2 e \left (2 c d+\sqrt {b^2-4 a c} e\right )+c \left (a \sqrt {b^2-4 a c} e^2-c d \left (\sqrt {b^2-4 a c} d-4 a e\right )\right )+b c \left (c d^2+e \left (2 \sqrt {b^2-4 a c} d-3 a e\right )\right )\right ) \tanh ^{-1}\left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {2 c d-\left (b-\sqrt {b^2-4 a c}\right ) e}}\right )}{c^{5/2} \sqrt {b^2-4 a c} \sqrt {2 c d-\left (b-\sqrt {b^2-4 a c}\right ) e}}-\frac {\sqrt {2} \left (b^3 e^2-b^2 e \left (2 c d-\sqrt {b^2-4 a c} e\right )+b c \left (c d^2-e \left (2 \sqrt {b^2-4 a c} d+3 a e\right )\right )-c \left (a \sqrt {b^2-4 a c} e^2-c d \left (\sqrt {b^2-4 a c} d+4 a e\right )\right )\right ) \tanh ^{-1}\left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {2 c d-\left (b+\sqrt {b^2-4 a c}\right ) e}}\right )}{c^{5/2} \sqrt {b^2-4 a c} \sqrt {2 c d-\left (b+\sqrt {b^2-4 a c}\right ) e}} \\ \end{align*}

Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 1.75 (sec) , antiderivative size = 493, normalized size of antiderivative = 1.09 \[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\frac {2 \sqrt {c} \sqrt {d+e x} (4 c d-3 b e+c e x)+\frac {3 \left (i b^3 e^2+b^2 e \left (-2 i c d+\sqrt {-b^2+4 a c} e\right )+i b c \left (c d^2+e \left (2 i \sqrt {-b^2+4 a c} d-3 a e\right )\right )+c \left (-a \sqrt {-b^2+4 a c} e^2+c d \left (\sqrt {-b^2+4 a c} d+4 i a e\right )\right )\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {-2 c d+b e-i \sqrt {-b^2+4 a c} e}}\right )}{\sqrt {-\frac {b^2}{2}+2 a c} \sqrt {-2 c d+\left (b-i \sqrt {-b^2+4 a c}\right ) e}}+\frac {3 \left (-i b^3 e^2+b^2 e \left (2 i c d+\sqrt {-b^2+4 a c} e\right )+b c \left (-i c d^2+e \left (-2 \sqrt {-b^2+4 a c} d+3 i a e\right )\right )+c \left (-a \sqrt {-b^2+4 a c} e^2+c d \left (\sqrt {-b^2+4 a c} d-4 i a e\right )\right )\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} \sqrt {d+e x}}{\sqrt {-2 c d+b e+i \sqrt {-b^2+4 a c} e}}\right )}{\sqrt {-\frac {b^2}{2}+2 a c} \sqrt {-2 c d+\left (b+i \sqrt {-b^2+4 a c}\right ) e}}}{3 c^{5/2}} \]

[In]

Integrate[(x*(d + e*x)^(3/2))/(a + b*x + c*x^2),x]

[Out]

(2*Sqrt[c]*Sqrt[d + e*x]*(4*c*d - 3*b*e + c*e*x) + (3*(I*b^3*e^2 + b^2*e*((-2*I)*c*d + Sqrt[-b^2 + 4*a*c]*e) +
 I*b*c*(c*d^2 + e*((2*I)*Sqrt[-b^2 + 4*a*c]*d - 3*a*e)) + c*(-(a*Sqrt[-b^2 + 4*a*c]*e^2) + c*d*(Sqrt[-b^2 + 4*
a*c]*d + (4*I)*a*e)))*ArcTan[(Sqrt[2]*Sqrt[c]*Sqrt[d + e*x])/Sqrt[-2*c*d + b*e - I*Sqrt[-b^2 + 4*a*c]*e]])/(Sq
rt[-1/2*b^2 + 2*a*c]*Sqrt[-2*c*d + (b - I*Sqrt[-b^2 + 4*a*c])*e]) + (3*((-I)*b^3*e^2 + b^2*e*((2*I)*c*d + Sqrt
[-b^2 + 4*a*c]*e) + b*c*((-I)*c*d^2 + e*(-2*Sqrt[-b^2 + 4*a*c]*d + (3*I)*a*e)) + c*(-(a*Sqrt[-b^2 + 4*a*c]*e^2
) + c*d*(Sqrt[-b^2 + 4*a*c]*d - (4*I)*a*e)))*ArcTan[(Sqrt[2]*Sqrt[c]*Sqrt[d + e*x])/Sqrt[-2*c*d + b*e + I*Sqrt
[-b^2 + 4*a*c]*e]])/(Sqrt[-1/2*b^2 + 2*a*c]*Sqrt[-2*c*d + (b + I*Sqrt[-b^2 + 4*a*c])*e]))/(3*c^(5/2))

Maple [A] (verified)

Time = 0.70 (sec) , antiderivative size = 473, normalized size of antiderivative = 1.04

method result size
pseudoelliptic \(-\frac {-\sqrt {2}\, \sqrt {\left (b e -2 c d +\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\right ) c}\, \left (\left (\left (a c -b^{2}\right ) e^{2}+2 b c d e -c^{2} d^{2}\right ) \sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}+\left (-3 a b c +b^{3}\right ) e^{3}+2 d \left (2 a \,c^{2}-b^{2} c \right ) e^{2}+b \,c^{2} d^{2} e \right ) \operatorname {arctanh}\left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (-b e +2 c d +\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\right ) c}}\right )+\left (\left (\left (\left (a c -b^{2}\right ) e^{2}+2 b c d e -c^{2} d^{2}\right ) \sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}+\left (3 a b c -b^{3}\right ) e^{3}+2 d \left (-2 a \,c^{2}+b^{2} c \right ) e^{2}-b \,c^{2} d^{2} e \right ) \sqrt {2}\, \arctan \left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (b e -2 c d +\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\right ) c}}\right )+2 \sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\, \sqrt {\left (b e -2 c d +\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\right ) c}\, \sqrt {e x +d}\, \left (\left (-\frac {c x}{3}+b \right ) e -\frac {4 c d}{3}\right )\right ) \sqrt {\left (-b e +2 c d +\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\right ) c}}{\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\, \sqrt {\left (b e -2 c d +\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\right ) c}\, \sqrt {\left (-b e +2 c d +\sqrt {-4 e^{2} \left (a c -\frac {b^{2}}{4}\right )}\right ) c}\, c^{2}}\) \(473\)
risch \(-\frac {2 \left (-x c e +3 b e -4 c d \right ) \sqrt {e x +d}}{3 c^{2}}-\frac {8 \left (-\frac {\left (-3 c \,e^{3} b a +4 a \,c^{2} d \,e^{2}+b^{3} e^{3}-2 b^{2} c d \,e^{2}+b \,c^{2} d^{2} e +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, a c \,e^{2}-\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b^{2} e^{2}+2 \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b c d e -\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, c^{2} d^{2}\right ) \sqrt {2}\, \operatorname {arctanh}\left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (-b e +2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}\right )}{8 c \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, \sqrt {\left (-b e +2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}+\frac {\left (3 c \,e^{3} b a -4 a \,c^{2} d \,e^{2}-b^{3} e^{3}+2 b^{2} c d \,e^{2}-b \,c^{2} d^{2} e +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, a c \,e^{2}-\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b^{2} e^{2}+2 \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b c d e -\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, c^{2} d^{2}\right ) \sqrt {2}\, \arctan \left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (b e -2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}\right )}{8 c \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, \sqrt {\left (b e -2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}\right )}{c}\) \(512\)
derivativedivides \(-\frac {2 \left (-\frac {c \left (e x +d \right )^{\frac {3}{2}}}{3}+b e \sqrt {e x +d}-c d \sqrt {e x +d}\right )}{c^{2}}+\frac {\frac {\left (-3 c \,e^{3} b a +4 a \,c^{2} d \,e^{2}+b^{3} e^{3}-2 b^{2} c d \,e^{2}+b \,c^{2} d^{2} e -\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, a c \,e^{2}+\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b^{2} e^{2}-2 \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b c d e +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, c^{2} d^{2}\right ) \sqrt {2}\, \arctan \left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (b e -2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}\right )}{c \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, \sqrt {\left (b e -2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}-\frac {\left (3 c \,e^{3} b a -4 a \,c^{2} d \,e^{2}-b^{3} e^{3}+2 b^{2} c d \,e^{2}-b \,c^{2} d^{2} e -\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, a c \,e^{2}+\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b^{2} e^{2}-2 \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b c d e +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, c^{2} d^{2}\right ) \sqrt {2}\, \operatorname {arctanh}\left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (-b e +2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}\right )}{c \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, \sqrt {\left (-b e +2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}}{c}\) \(521\)
default \(-\frac {2 \left (-\frac {c \left (e x +d \right )^{\frac {3}{2}}}{3}+b e \sqrt {e x +d}-c d \sqrt {e x +d}\right )}{c^{2}}+\frac {\frac {\left (-3 c \,e^{3} b a +4 a \,c^{2} d \,e^{2}+b^{3} e^{3}-2 b^{2} c d \,e^{2}+b \,c^{2} d^{2} e -\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, a c \,e^{2}+\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b^{2} e^{2}-2 \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b c d e +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, c^{2} d^{2}\right ) \sqrt {2}\, \arctan \left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (b e -2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}\right )}{c \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, \sqrt {\left (b e -2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}-\frac {\left (3 c \,e^{3} b a -4 a \,c^{2} d \,e^{2}-b^{3} e^{3}+2 b^{2} c d \,e^{2}-b \,c^{2} d^{2} e -\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, a c \,e^{2}+\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b^{2} e^{2}-2 \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, b c d e +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, c^{2} d^{2}\right ) \sqrt {2}\, \operatorname {arctanh}\left (\frac {c \sqrt {e x +d}\, \sqrt {2}}{\sqrt {\left (-b e +2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}\right )}{c \sqrt {-e^{2} \left (4 a c -b^{2}\right )}\, \sqrt {\left (-b e +2 c d +\sqrt {-e^{2} \left (4 a c -b^{2}\right )}\right ) c}}}{c}\) \(521\)

[In]

int(x*(e*x+d)^(3/2)/(c*x^2+b*x+a),x,method=_RETURNVERBOSE)

[Out]

-(-2^(1/2)*((b*e-2*c*d+(-4*e^2*(a*c-1/4*b^2))^(1/2))*c)^(1/2)*(((a*c-b^2)*e^2+2*b*c*d*e-c^2*d^2)*(-4*e^2*(a*c-
1/4*b^2))^(1/2)+(-3*a*b*c+b^3)*e^3+2*d*(2*a*c^2-b^2*c)*e^2+b*c^2*d^2*e)*arctanh(c*(e*x+d)^(1/2)*2^(1/2)/((-b*e
+2*c*d+(-4*e^2*(a*c-1/4*b^2))^(1/2))*c)^(1/2))+((((a*c-b^2)*e^2+2*b*c*d*e-c^2*d^2)*(-4*e^2*(a*c-1/4*b^2))^(1/2
)+(3*a*b*c-b^3)*e^3+2*d*(-2*a*c^2+b^2*c)*e^2-b*c^2*d^2*e)*2^(1/2)*arctan(c*(e*x+d)^(1/2)*2^(1/2)/((b*e-2*c*d+(
-4*e^2*(a*c-1/4*b^2))^(1/2))*c)^(1/2))+2*(-4*e^2*(a*c-1/4*b^2))^(1/2)*((b*e-2*c*d+(-4*e^2*(a*c-1/4*b^2))^(1/2)
)*c)^(1/2)*(e*x+d)^(1/2)*((-1/3*c*x+b)*e-4/3*c*d))*((-b*e+2*c*d+(-4*e^2*(a*c-1/4*b^2))^(1/2))*c)^(1/2))/(-4*e^
2*(a*c-1/4*b^2))^(1/2)/((b*e-2*c*d+(-4*e^2*(a*c-1/4*b^2))^(1/2))*c)^(1/2)/((-b*e+2*c*d+(-4*e^2*(a*c-1/4*b^2))^
(1/2))*c)^(1/2)/c^2

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 5572 vs. \(2 (392) = 784\).

Time = 1.65 (sec) , antiderivative size = 5572, normalized size of antiderivative = 12.30 \[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\text {Too large to display} \]

[In]

integrate(x*(e*x+d)^(3/2)/(c*x^2+b*x+a),x, algorithm="fricas")

[Out]

Too large to include

Sympy [F(-1)]

Timed out. \[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\text {Timed out} \]

[In]

integrate(x*(e*x+d)**(3/2)/(c*x**2+b*x+a),x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\int { \frac {{\left (e x + d\right )}^{\frac {3}{2}} x}{c x^{2} + b x + a} \,d x } \]

[In]

integrate(x*(e*x+d)^(3/2)/(c*x^2+b*x+a),x, algorithm="maxima")

[Out]

integrate((e*x + d)^(3/2)*x/(c*x^2 + b*x + a), x)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 986 vs. \(2 (392) = 784\).

Time = 0.36 (sec) , antiderivative size = 986, normalized size of antiderivative = 2.18 \[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\frac {2 \, {\left ({\left (e x + d\right )}^{\frac {3}{2}} c^{2} + 3 \, \sqrt {e x + d} c^{2} d - 3 \, \sqrt {e x + d} b c e\right )}}{3 \, c^{3}} + \frac {{\left (\sqrt {-4 \, c^{2} d + 2 \, {\left (b c - \sqrt {b^{2} - 4 \, a c} c\right )} e} {\left ({\left (b^{2} c^{2} - 4 \, a c^{3}\right )} d^{2} - 2 \, {\left (b^{3} c - 4 \, a b c^{2}\right )} d e + {\left (b^{4} - 5 \, a b^{2} c + 4 \, a^{2} c^{2}\right )} e^{2}\right )} c^{2} e^{2} - 2 \, {\left (\sqrt {b^{2} - 4 \, a c} c^{4} d^{3} - 2 \, \sqrt {b^{2} - 4 \, a c} b c^{3} d^{2} e - \sqrt {b^{2} - 4 \, a c} a b c^{2} e^{3} + {\left (b^{2} c^{2} + a c^{3}\right )} \sqrt {b^{2} - 4 \, a c} d e^{2}\right )} \sqrt {-4 \, c^{2} d + 2 \, {\left (b c - \sqrt {b^{2} - 4 \, a c} c\right )} e} {\left | c \right |} {\left | e \right |} + {\left (2 \, b c^{5} d^{3} e - {\left (5 \, b^{2} c^{4} - 8 \, a c^{5}\right )} d^{2} e^{2} + 2 \, {\left (2 \, b^{3} c^{3} - 5 \, a b c^{4}\right )} d e^{3} - {\left (b^{4} c^{2} - 3 \, a b^{2} c^{3}\right )} e^{4}\right )} \sqrt {-4 \, c^{2} d + 2 \, {\left (b c - \sqrt {b^{2} - 4 \, a c} c\right )} e}\right )} \arctan \left (\frac {2 \, \sqrt {\frac {1}{2}} \sqrt {e x + d}}{\sqrt {-\frac {2 \, c^{4} d - b c^{3} e + \sqrt {-4 \, {\left (c^{4} d^{2} - b c^{3} d e + a c^{3} e^{2}\right )} c^{4} + {\left (2 \, c^{4} d - b c^{3} e\right )}^{2}}}{c^{4}}}}\right )}{4 \, {\left (\sqrt {b^{2} - 4 \, a c} c^{5} d^{2} - \sqrt {b^{2} - 4 \, a c} b c^{4} d e + \sqrt {b^{2} - 4 \, a c} a c^{4} e^{2}\right )} c^{2} {\left | e \right |}} - \frac {{\left (\sqrt {-4 \, c^{2} d + 2 \, {\left (b c + \sqrt {b^{2} - 4 \, a c} c\right )} e} {\left ({\left (b^{2} c^{2} - 4 \, a c^{3}\right )} d^{2} - 2 \, {\left (b^{3} c - 4 \, a b c^{2}\right )} d e + {\left (b^{4} - 5 \, a b^{2} c + 4 \, a^{2} c^{2}\right )} e^{2}\right )} c^{2} e^{2} + 2 \, {\left (\sqrt {b^{2} - 4 \, a c} c^{4} d^{3} - 2 \, \sqrt {b^{2} - 4 \, a c} b c^{3} d^{2} e - \sqrt {b^{2} - 4 \, a c} a b c^{2} e^{3} + {\left (b^{2} c^{2} + a c^{3}\right )} \sqrt {b^{2} - 4 \, a c} d e^{2}\right )} \sqrt {-4 \, c^{2} d + 2 \, {\left (b c + \sqrt {b^{2} - 4 \, a c} c\right )} e} {\left | c \right |} {\left | e \right |} + {\left (2 \, b c^{5} d^{3} e - {\left (5 \, b^{2} c^{4} - 8 \, a c^{5}\right )} d^{2} e^{2} + 2 \, {\left (2 \, b^{3} c^{3} - 5 \, a b c^{4}\right )} d e^{3} - {\left (b^{4} c^{2} - 3 \, a b^{2} c^{3}\right )} e^{4}\right )} \sqrt {-4 \, c^{2} d + 2 \, {\left (b c + \sqrt {b^{2} - 4 \, a c} c\right )} e}\right )} \arctan \left (\frac {2 \, \sqrt {\frac {1}{2}} \sqrt {e x + d}}{\sqrt {-\frac {2 \, c^{4} d - b c^{3} e - \sqrt {-4 \, {\left (c^{4} d^{2} - b c^{3} d e + a c^{3} e^{2}\right )} c^{4} + {\left (2 \, c^{4} d - b c^{3} e\right )}^{2}}}{c^{4}}}}\right )}{4 \, {\left (\sqrt {b^{2} - 4 \, a c} c^{5} d^{2} - \sqrt {b^{2} - 4 \, a c} b c^{4} d e + \sqrt {b^{2} - 4 \, a c} a c^{4} e^{2}\right )} c^{2} {\left | e \right |}} \]

[In]

integrate(x*(e*x+d)^(3/2)/(c*x^2+b*x+a),x, algorithm="giac")

[Out]

2/3*((e*x + d)^(3/2)*c^2 + 3*sqrt(e*x + d)*c^2*d - 3*sqrt(e*x + d)*b*c*e)/c^3 + 1/4*(sqrt(-4*c^2*d + 2*(b*c -
sqrt(b^2 - 4*a*c)*c)*e)*((b^2*c^2 - 4*a*c^3)*d^2 - 2*(b^3*c - 4*a*b*c^2)*d*e + (b^4 - 5*a*b^2*c + 4*a^2*c^2)*e
^2)*c^2*e^2 - 2*(sqrt(b^2 - 4*a*c)*c^4*d^3 - 2*sqrt(b^2 - 4*a*c)*b*c^3*d^2*e - sqrt(b^2 - 4*a*c)*a*b*c^2*e^3 +
 (b^2*c^2 + a*c^3)*sqrt(b^2 - 4*a*c)*d*e^2)*sqrt(-4*c^2*d + 2*(b*c - sqrt(b^2 - 4*a*c)*c)*e)*abs(c)*abs(e) + (
2*b*c^5*d^3*e - (5*b^2*c^4 - 8*a*c^5)*d^2*e^2 + 2*(2*b^3*c^3 - 5*a*b*c^4)*d*e^3 - (b^4*c^2 - 3*a*b^2*c^3)*e^4)
*sqrt(-4*c^2*d + 2*(b*c - sqrt(b^2 - 4*a*c)*c)*e))*arctan(2*sqrt(1/2)*sqrt(e*x + d)/sqrt(-(2*c^4*d - b*c^3*e +
 sqrt(-4*(c^4*d^2 - b*c^3*d*e + a*c^3*e^2)*c^4 + (2*c^4*d - b*c^3*e)^2))/c^4))/((sqrt(b^2 - 4*a*c)*c^5*d^2 - s
qrt(b^2 - 4*a*c)*b*c^4*d*e + sqrt(b^2 - 4*a*c)*a*c^4*e^2)*c^2*abs(e)) - 1/4*(sqrt(-4*c^2*d + 2*(b*c + sqrt(b^2
 - 4*a*c)*c)*e)*((b^2*c^2 - 4*a*c^3)*d^2 - 2*(b^3*c - 4*a*b*c^2)*d*e + (b^4 - 5*a*b^2*c + 4*a^2*c^2)*e^2)*c^2*
e^2 + 2*(sqrt(b^2 - 4*a*c)*c^4*d^3 - 2*sqrt(b^2 - 4*a*c)*b*c^3*d^2*e - sqrt(b^2 - 4*a*c)*a*b*c^2*e^3 + (b^2*c^
2 + a*c^3)*sqrt(b^2 - 4*a*c)*d*e^2)*sqrt(-4*c^2*d + 2*(b*c + sqrt(b^2 - 4*a*c)*c)*e)*abs(c)*abs(e) + (2*b*c^5*
d^3*e - (5*b^2*c^4 - 8*a*c^5)*d^2*e^2 + 2*(2*b^3*c^3 - 5*a*b*c^4)*d*e^3 - (b^4*c^2 - 3*a*b^2*c^3)*e^4)*sqrt(-4
*c^2*d + 2*(b*c + sqrt(b^2 - 4*a*c)*c)*e))*arctan(2*sqrt(1/2)*sqrt(e*x + d)/sqrt(-(2*c^4*d - b*c^3*e - sqrt(-4
*(c^4*d^2 - b*c^3*d*e + a*c^3*e^2)*c^4 + (2*c^4*d - b*c^3*e)^2))/c^4))/((sqrt(b^2 - 4*a*c)*c^5*d^2 - sqrt(b^2
- 4*a*c)*b*c^4*d*e + sqrt(b^2 - 4*a*c)*a*c^4*e^2)*c^2*abs(e))

Mupad [B] (verification not implemented)

Time = 13.80 (sec) , antiderivative size = 13841, normalized size of antiderivative = 30.55 \[ \int \frac {x (d+e x)^{3/2}}{a+b x+c x^2} \, dx=\text {Too large to display} \]

[In]

int((x*(d + e*x)^(3/2))/(a + b*x + c*x^2),x)

[Out]

(2*(d + e*x)^(3/2))/(3*c) - ((2*d)/c + (2*(b*e - 2*c*d))/c^2)*(d + e*x)^(1/2) - atan(((((8*(a*b^3*c^3*e^5 - 4*
a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^3*e^2 + 2*b^3*c^4*d^2*e^3 - 8*a*
b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 - (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^
3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*
c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^
6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^
2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 +
 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5
 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c^7*d*e^2))/c^3)*(-(b^7*e^3 - 8*a
^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(
-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2
)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*
a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^
3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3
)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) - (8*(d + e*x)^(1/2)*(b^6*e^6 - 2*a^3*c^3*e^6 - 2*a*c
^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*c^3*d^3*e^3 + 6*b^4*c^2*d^2*e^4
- 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20*a^2*b*c^3*d*e^5 - 24*a*b^2*c^3*
d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 -
 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e
^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(
1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) -
 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*
a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*1i - (((8*(a*b^3*c^3*e
^5 - 4*a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^3*e^2 + 2*b^3*c^4*d^2*e^3
 - 8*a*b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 + (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 +
 b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) +
24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3
 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b
*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*
d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 +
b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c^7*d*e^2))/c^3)*(-(b^7*e^
3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^
3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*
c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*
e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c
- b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c -
 b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) + (8*(d + e*x)^(1/2)*(b^6*e^6 - 2*a^3*c^3*e^6
- 2*a*c^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*c^3*d^3*e^3 + 6*b^4*c^2*d
^2*e^4 - 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20*a^2*b*c^3*d*e^5 - 24*a*b
^2*c^3*d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^
4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^
3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^
2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^
(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/
2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*1i)/((16*(a^4*c
*e^8 - a^3*b^2*e^8 - a*b^4*d^2*e^6 + 2*a^2*b^3*d*e^7 - a*c^4*d^6*e^2 - a^2*c^3*d^4*e^4 + a^3*c^2*d^2*e^6 + 4*a
*b*c^3*d^5*e^3 + 4*a*b^3*c*d^3*e^5 - 6*a*b^2*c^2*d^4*e^4 + 4*a^2*b*c^2*d^3*e^5 - 5*a^2*b^2*c*d^2*e^6))/c^3 + (
((8*(a*b^3*c^3*e^5 - 4*a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^3*e^2 + 2
*b^3*c^4*d^2*e^3 - 8*a*b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 - (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^2*c^5*d^3
 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c -
b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2)
 - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*
d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) -
 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(
2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c^7*d*e^2)
)/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*
b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2
*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 2
1*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c
*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*
d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) - (8*(d + e*x)^(1/2)*(b^6*e^6
- 2*a^3*c^3*e^6 - 2*a*c^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*c^3*d^3*e
^3 + 6*b^4*c^2*d^2*e^4 - 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20*a^2*b*c^
3*d*e^5 - 24*a*b^2*c^3*d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1
/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d
^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e
^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(
4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a
*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)
 + (((8*(a*b^3*c^3*e^5 - 4*a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^3*e^2
 + 2*b^3*c^4*d^2*e^3 - 8*a*b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 + (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^2*c^5
*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*
c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(
1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*
c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/
2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2
))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c^7*d*
e^2))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*
a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 +
 a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2)
 - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b
^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*
c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) + (8*(d + e*x)^(1/2)*(b^6*
e^6 - 2*a^3*c^3*e^6 - 2*a*c^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*c^3*d
^3*e^3 + 6*b^4*c^2*d^2*e^4 - 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20*a^2*
b*c^3*d*e^5 - 24*a*b^2*c^3*d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3
)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c
^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2
*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e
*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-
(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(
1/2)))*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 + b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*
b*c^3*e^3 - b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 + a^2
*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 - 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 2
1*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e - 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 3*b^3*c
*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 + 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b*c^2*
d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*2i - atan(((((8*(a*b^3*c^3*e^5
 - 4*a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^3*e^2 + 2*b^3*c^4*d^2*e^3 -
 8*a*b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 - (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b
^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24
*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 -
 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c
^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*
e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^
4*c^5 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c^7*d*e^2))/c^3)*(-(b^7*e^3
- 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*
d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c
- b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e
+ 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c -
b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b
^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) - (8*(d + e*x)^(1/2)*(b^6*e^6 - 2*a^3*c^3*e^6 -
2*a*c^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*c^3*d^3*e^3 + 6*b^4*c^2*d^2
*e^4 - 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20*a^2*b*c^3*d*e^5 - 24*a*b^2
*c^3*d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*
d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*
c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)
^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1
/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2)
 - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*1i - (((8*(a*b^3*
c^3*e^5 - 4*a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^3*e^2 + 2*b^3*c^4*d^
2*e^3 - 8*a*b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 + (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*
d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/
2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*
c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*
a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2
*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c
^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c^7*d*e^2))/c^3)*(-(b
^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 +
 b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-
(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3
*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4
*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*
a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) + (8*(d + e*x)^(1/2)*(b^6*e^6 - 2*a^3*c^3
*e^6 - 2*a*c^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*c^3*d^3*e^3 + 6*b^4*
c^2*d^2*e^4 - 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20*a^2*b*c^3*d*e^5 - 2
4*a*b^2*c^3*d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b
^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a
^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c
 - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2
)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3
)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*1i)/((16*(
a^4*c*e^8 - a^3*b^2*e^8 - a*b^4*d^2*e^6 + 2*a^2*b^3*d*e^7 - a*c^4*d^6*e^2 - a^2*c^3*d^4*e^4 + a^3*c^2*d^2*e^6
+ 4*a*b*c^3*d^5*e^3 + 4*a*b^3*c*d^3*e^5 - 6*a*b^2*c^2*d^4*e^4 + 4*a^2*b*c^2*d^3*e^5 - 5*a^2*b^2*c*d^2*e^6))/c^
3 + (((8*(a*b^3*c^3*e^5 - 4*a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^3*e^
2 + 2*b^3*c^4*d^2*e^3 - 8*a*b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 - (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^2*c^
5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a
*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^
(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4
*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1
/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/
2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c^7*d
*e^2))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20
*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3
- a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2
) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*
b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b
*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) - (8*(d + e*x)^(1/2)*(b^6
*e^6 - 2*a^3*c^3*e^6 - 2*a*c^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*c^3*
d^3*e^3 + 6*b^4*c^2*d^2*e^4 - 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20*a^2
*b*c^3*d*e^5 - 24*a*b^2*c^3*d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^
3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*
c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^
2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*
e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(
-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^
(1/2) + (((8*(a*b^3*c^3*e^5 - 4*a^2*b*c^4*e^5 + 4*a*c^6*d^3*e^2 + 4*a^2*c^5*d*e^4 - b^4*c^3*d*e^4 - b^2*c^5*d^
3*e^2 + 2*b^3*c^4*d^2*e^3 - 8*a*b*c^5*d^2*e^3 + 3*a*b^2*c^4*d*e^4))/c^3 + (8*(d + e*x)^(1/2)*(-(b^7*e^3 - 8*a^
2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-
(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)
^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a
*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3
)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)
^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*(b^3*c^5*e^3 - 2*b^2*c^6*d*e^2 - 4*a*b*c^6*e^3 + 8*a*c
^7*d*e^2))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3
- 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*
e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^
(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2)
+ 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6
*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2) + (8*(d + e*x)^(1/2)*
(b^6*e^6 - 2*a^3*c^3*e^6 - 2*a*c^5*d^4*e^2 + 9*a^2*b^2*c^2*e^6 + 12*a^2*c^4*d^2*e^4 + b^2*c^4*d^4*e^2 - 4*b^3*
c^3*d^3*e^3 + 6*b^4*c^2*d^2*e^4 - 6*a*b^4*c*e^6 - 4*b^5*c*d*e^5 + 12*a*b*c^4*d^3*e^3 + 20*a*b^3*c^2*d*e^5 - 20
*a^2*b*c^3*d*e^5 - 24*a*b^2*c^3*d^2*e^4))/c^3)*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b
^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*
b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3 - a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*
a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*
d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2
*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6
)))^(1/2)))*(-(b^7*e^3 - 8*a^2*c^5*d^3 - b^4*c^3*d^3 - b^4*e^3*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^2*c^4*d^3 - 20
*a^3*b*c^3*e^3 + b*c^3*d^3*(-(4*a*c - b^2)^3)^(1/2) + 24*a^3*c^4*d*e^2 + 3*b^5*c^2*d^2*e + 25*a^2*b^3*c^2*e^3
- a^2*c^2*e^3*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*e^3 - 3*b^6*c*d*e^2 + 3*a*b^2*c*e^3*(-(4*a*c - b^2)^3)^(1/2
) - 21*a*b^3*c^3*d^2*e + 24*a*b^4*c^2*d*e^2 + 36*a^2*b*c^4*d^2*e + 3*a*c^3*d^2*e*(-(4*a*c - b^2)^3)^(1/2) + 3*
b^3*c*d*e^2*(-(4*a*c - b^2)^3)^(1/2) - 54*a^2*b^2*c^3*d*e^2 - 3*b^2*c^2*d^2*e*(-(4*a*c - b^2)^3)^(1/2) - 6*a*b
*c^2*d*e^2*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^7 + b^4*c^5 - 8*a*b^2*c^6)))^(1/2)*2i